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Abstract 

Diagnosis  and  prognosis  are  necessary  tasks  for  system  re¬ 
configuration  and  fault-adaptive  control  in  complex  systems. 
Diagnosis  consists  of  detection,  isolation  and  identification  of 
faults,  while  prognosis  consists  of  prediction  of  the  remain¬ 
ing  useful  life  of  systems.  This  paper  presents  a  novel  inte¬ 
grated  framework  for  model-based  distributed  diagnosis  and 
prognosis,  where  system  decomposition  is  used  to  enable  the 
diagnosis  and  prognosis  tasks  to  be  performed  in  a  distributed 
way.  We  show  how  different  submodels  can  be  automati¬ 
cally  constructed  to  solve  the  local  diagnosis  and  prognosis 
problems.  We  illustrate  our  approach  using  a  simulated  four- 
wheeled  rover  for  different  fault  scenarios.  Our  experiments 
show  that  our  approach  correctly  performs  distributed  fault 
diagnosis  and  prognosis  in  an  efficient  and  robust  manner. 

1.  Introduction 

Systems  health  monitoring  is  essential  to  guaranteeing  the 
safe,  efficient,  and  reliable  operation  of  engineering  sys¬ 
tems.  Integrated  systems  health  management  methodologies 
include  fault  diagnosis  and  prognosis  mechanisms,  where  di¬ 
agnosis  involves  detecting  when  a  fault  has  occurred,  isolat¬ 
ing  the  true  fault,  and  identifying  the  true  damage  to  the  sys¬ 
tem;  and  prognosis  involves  predicting  how  much  useful  life 
remains  in  the  different  components,  subsystems,  or  systems 
given  the  diagnosed  fault  conditions.  The  information  on  the 
fault  size  and  its  expected  impact  on  system  life  can  be  used 
to  initiate  recovery  and  reconfiguration  actions  that  mitigate 
the  fault  or  extend  system  life. 
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A  large  body  of  research  exists  for  both  model-based  di¬ 
agnosis  (Gertler,  1998;  Blanke  et  al.,  2006)  and  prognosis 
methods  (Luo  et  al.,  2008;  Saha  &  Goebel,  2009;  Orchard 
&  Vachtsevanos,  2009).  However,  the  integration  of  diag¬ 
nosis  and  prognosis  algorithms  is  seldom  studied.  In  fact, 
many  diagnosis  methodologies  leave  out  the  fault  identifica¬ 
tion  step,  which  is  necessary  to  perform  a  prediction  from  the 
current  system  state.  Recently,  we  presented  an  integrated 
model-based  framework  for  diagnosis  and  prognosis  of  com¬ 
plex  systems,  in  which  we  made  use  of  a  common  modeling 
framework  for  modeling  both  the  nominal  and  faulty  system 
behavior  (Roychoudhury  &  Daigle,  201 1). 

In  (Roychoudhury  &  Daigle,  2011),  the  nominal  system  be¬ 
havior  is  estimated  using  an  observer  built  with  the  nominal 
model.  Faults  are  detected  when  a  statistically  significant  de¬ 
viation  between  the  nominal  estimates  and  the  observed  mea¬ 
surements  is  observed  (Biswas  et  al.,  2003).  Fault  isolation 
compares  the  observed  measurement  deviations  against  pre¬ 
dictions  of  how  the  measurements  would  deviate  for  each 
possible  fault  (Mosterman  &  Biswas,  1999).  Fault  identifi¬ 
cation  performs  joint  state-parameter  estimation  using  multi¬ 
ple  observers,  where,  for  each  fault,  the  faulty  system  model 
is  constmcted  as  the  nominal  model  integrated  with  a  hy¬ 
pothesized  fault  model  (Roychoudhury,  2009).  The  prog¬ 
nosis  module  uses,  for  each  fault  hypothesis,  a  prediction 
model  based  on  its  faulty  system  model  and  the  identified 
fault  parameters,  to  predict  the  remaining  useful  life  of  the 
system  (Daigle,  Saha,  &  Goebel,  2012).  However,  this  inte¬ 
grated  solution  performs  the  diagnosis  and  prognosis  task  in  a 
centralized  fashion,  which  is  prone  to  single  points  of  failure, 
and  does  not  scale  well  as  the  size  of  the  system  increases. 

To  overcome  such  problems,  in  this  work,  we  leverage  re¬ 
cent  results  for  distributed  diagnosis  (Bregon  et  al.,  2011) 
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and  distributed  prognosis  (Daigle,  Bregon,  &  Roychoud- 
hury,  2012),  which  make  use  of  structural  model  decompo¬ 
sition  techniques,  to  provide  a  systematic  approach  to  dis¬ 
tributing  the  different  diagnosis  and  prognosis  steps  presented 
in  (Roychoudhury  &  Daigle,  201 1). 

Distributed  diagnosis  is  achieved  by  designing  local  dis¬ 
tributed  subsystems  based  on  global  diagnosability  analysis 
of  the  system,  thus  computing  globally  correct  distributed 
diagnosis  results  without  the  use  of  a  centralized  coordina¬ 
tor  (Bregon  et  al.,  2011).  These  local  distributed  subsystems 
are  then  used  to  construct  local  event-based  distributed  diag¬ 
noses  for  distributed  fault  isolation.  Distributed  fault  iden¬ 
tification  is  achieved  by  developing  independent  local  state- 
parameter  estimators  for  each  hypothesized  fault.  Regarding 
distributed  prediction,  in  (Daigle,  Bregon,  &  Roychoudhury, 
2012)  we  developed  an  architecture  that  enables  a  large  prog¬ 
nosis  problem  to  be  decomposed  into  several  independent  lo¬ 
cal  subproblems  from  which  local  results  can  be  merged  into 
a  global  result. 

The  main  contribution  of  this  paper  is  an  integrated  frame¬ 
work  for  distributed  model-based  diagnosis  and  prognosis  of 
single  faults  based  on  structural  model  decomposition.  The 
proposed  framework  scales  well  and  the  resulting  subprob¬ 
lems  are  typically  small  and  easy  to  solve,  resulting  in  an  ef¬ 
ficient  and  scalable  distributed  solution  to  the  combined  diag¬ 
nosis  and  prognosis  problem.  We  perform  a  number  of  exper¬ 
iments  on  a  simulated  four-wheeled  rover  testbed  (Balaban  et 
al.,  201 1)  to  demonstrate  and  evaluate  our  approach. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2  pro¬ 
vides  the  problem  formulation  for  our  diagnosis  and  progno¬ 
sis  framework.  Section  3  describes  the  distributed  architec¬ 
ture  and  Section  4  briefly  introduces  its  different  components. 
Section  5  presents  the  case  study  and  experimental  results. 
Finally,  Section  6  concludes  the  paper. 

2.  Problem  Formulation 

The  nominal  system  model  is  represented  as  follows: 

*(<)  =  f(t,x(t),0(t),u(t),v(t)), 
y(t)  =  h(f,x(f),0(f),u(f),n(f)), 

where  x(f)  G  R”*  is  the  state  vector,  0(f)  G  R”8  is  the 
parameter  vector,  u(f)  G  R"“  is  the  input  vector,  v(t)  G  R”7' 
is  the  process  noise  vector,  f  is  the  state  equation,  y(t)  G  R"" 
is  the  output  vector,  n(f)  G  R"71  is  the  measurement  noise 
vector,  and  h  is  the  output  equation. 1 

Faults  in  the  system  are  represented  as  changes  in  the  above 
nominal  system  model.  In  this  work,  we  only  consider  sin¬ 
gle  faults  occurring  as  changes  in  system  parameters,  0(f). 
We  denote  a  fault,  /  G  F,  as  a  tuple,  ( 9,gy ),  where,  0  G  0 

'Here,  we  use  bold  typeface  to  denote  vectors,  and  use  na  to  denote  the 
length  of  a  vector  a. 


is  the  fault  parameter ,  and  g /  denotes  the  fault  progression 
function,  which  models  the  way  fault  /  is  manifested  in  pa¬ 
rameter  9,  i.e., 

0(f)  =  0/(f,x/(f),0/(f),u(f),m/(f)), 

where  x/(f)  =  [x(f),  0(f)]T,  0/(f)  =  [0(f)\{0(f)}, 

4>f(t)]T ,  4>f(t)  G  M.n*f  is  a  vector  of  fault  progression  pa¬ 
rameters,  and  m/(f)  G  R nmf  is  a  process  noise  vector  asso¬ 
ciated  with  the  fault  progression  function. 

To  develop  our  integrated  diagnosis  and  prognosis  frame¬ 
work,  the  faulty  system  model  for  fault  /  =  (0,5/)  is  con¬ 
structed  from  the  nominal  system  model  by  including  the  pa¬ 
rameter  as  a  state  and  augmenting  the  state  equation  by  in¬ 
cluding  the  fault  progression  function,  i.e., 

X7W  =  f/(7x/(f),0/(f),u(f),v(f)), 
y(f)  =  h(f,  x(f),  0(f),  u(f),  n(f)), 

where, 

f  (.1  =  f  f(f,x(f),0(f),u(f),v(f))  i  =  r  x(f) 

/u  L  J  [  0(f) 

The  goal  of  diagnosis  is  to:  (i)  detect  a  change  in  some  0  G  0; 
(ii)  isolate,  under  the  single  fault  assumption,  the  true  fault 
/  G  F,  i.e.,  both  the  parameter  0  that  has  changed,  and  its 
fault  progression  function  gj;  and  (Hi)  identify  (i.e.  estimate) 
the  fault  by  computing  p(x/ (f) ,  0/(f)|y(O:f)),  where  y(0:f) 
denotes  all  measurements  observed  up  to  time  f. 

The  goal  of  prognosis  is  to  determine  the  end  of  (use¬ 
ful)  life  (EOL)  of  a  system,  and/or  its  remaining  useful 
life  (RUL).  For  a  given  fault,  /,  using  the  fault  estimate, 
p(x/(f),  0/(f)|y(O  :f)),  a  probability  distribution  of  EOL, 
p(EOL/(fP) |y(0 :  fP)),  and/or  RUL,  p(RUL/(fP)  |y(0 :  tP)) 
is  computed  at  a  given  time  point  tp  (Daigle,  Saha,  & 
Goebel,  2012).  Since  the  prognosis  problem  is  stochas¬ 
tic,  EOL/RUL  are  random  variables  and  we  represent  them 
by  probability  distributions.  The  acceptable  behavior  of 
the  system  is  expressed  through  a  set  of  nc  constraints, 
CEOLf  =  {cj}”=  i,  where  Cj  :  R"1/  x  R"8r  x  R"u  — >  B 
maps  a  given  point  in  the  joint  state-parameter  space  given 
the  current  inputs,  (x/(f),  0/(f),  u(f)),  to  the  Boolean  do¬ 
main  B  =  [0,1],  where  Ci(x/(f),  0/(f),  u(f))  =  1  if  the 
constraint  is  satisfied  (Daigle,  Saha,  &  Goebel,  2012).  If 
c,;(x/(f),  0/(f),  u(f))  =  0,  then  the  constraint  is  not  satis¬ 
fied,  and  the  behavior  of  the  system  is  deemed  to  be  unac¬ 
ceptable.  These  individual  constraints  are  combined  into  a 
single  threshold  function  Teol/  :  Rn:E/  x  R"8/  x  R"u  — >  B, 
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defined  as 

TEOLr(x/(f),0/(f),u(f))  = 

1,  0  G  {Ci(x/(f),  0/(f),  u^))}^ 

0,  otherwise. 

So,  EOL f  may  be  defined  as 

EOL  f(tP)  = 

inf{f  €R:t>tP  and  TEOL/  (x/(f),  0/(f),  u(i))  =  1}, 

i.e.,  EOL  is  the  earliest  time  point  at  which  the  threshold  is 
reached.  RUL  is  expressed  given  EOL  as 

RUL/(fP)  =  EOL f(tp)  -  tP. 

3.  Distributed  Architecture 

For  a  large  system,  both  the  diagnosis  and  prognosis  prob¬ 
lems  are  correspondingly  large.  A  centralized  approach  does 
not  scale  well,  can  be  computationally  expensive,  and  prone 
to  single  points  of  failure.  Therefore,  we  propose  to  decom¬ 
pose  the  global  integrated  diagnosis  and  prognosis  problem 
into  independent  local  subproblems.  In  this  work,  we  build 
on  the  ideas  from  structural  model  decomposition  (Blanke  et 
ah,  2006;  Pulido  &  Alonso-Gonzalez,  2004)  to  compute  lo¬ 
cal  independent  subproblems,  which  may  be  solved  in  paral¬ 
lel,  thus  providing  scalability  and  efficiency.  Model  decom¬ 
position  is  not  a  new  concept,  and  several  approaches  have 
been  developed  for  purposes  of  system  identification,  estima¬ 
tion,  learning,  and  diagnosis  (Staroswiecki  &  Declerck,  1989; 
Pulido  &  Alonso-Gonzalez,  2004;  Williams  &  Millar,  1998). 
Structural  model  decomposition  allows  decomposing  a  global 
model  into  a  set  of  local  submodels  for  which  local  diagnosis 
and  prognosis  problems  can  be  directly  defined.  The  global 
model  of  the  system,  denoted  as  M,  is  defined  as  follows. 

Definition  1  (Model).  The  model  of  a  system.  Ad,  is  a  tuple 
M  =  (X,  Q,U,Y,C ),  where  X  is  the  set  of  state  variables 
of  x,  0  is  the  set  of  unknown  parameters  of  0,  U  is  the  set  of 
input  variables  of  u,  Y  is  the  set  of  output  variables  of  y,  and 
C  is  the  set  of  model  constraints  of  f,  h,  and  EOL  constraints 
of  CpOLp 

The  basic  idea  of  the  model  decomposition  problem  is  to  de¬ 
compose  the  global  system  model  into  a  set  of  submodels  sat¬ 
isfying  given  constraints,  such  that  each  submodel  contains 
sufficient  analytical  redundancy  to  generate  fault  hypotheses 
from  observed  measurement  deviations.  A  submodel  in  our 
framework  is  defined  as  follows. 

Definition  2  (Submodel).  A  submodel  AT;  of  a  system  model 
M  =  (X,  0,  U,  Y,  C )  is  a  tuple  Mi  =  (X,,  0j}  Uu  Yit  Q), 
where  X*  C  X,  0.;  C  0,  Ui  C  X  U  U  U  Y,  and  Yt  C  Y  are 

the  state,  parameter,  input,  and  output  variables,  respectively. 


and  Ci  C  C  are  the  submodel  constraints.2 

Next,  we  discuss  the  fundamental  ideas  of  our  model  decom¬ 
position  approach  and  show  the  constraints  needed  to  obtain 
the  different  submodels  for  distributed  diagnosis  and  progno¬ 
sis.  Then,  we  propose  our  integrated  approach. 

3.1.  Model  Decomposition  for  Distributed  Diagnosis  and 
Prognosis 

Model  decomposition  in  our  framework  can  be  accomplished 
by  using  some  variables  (either  measured  variables  or  vari¬ 
ables  for  which  the  values  are  known)  as  local  inputs,  Ui ,  such 
that  each  one  of  the  submodels  satisfies  a  set  of  constraints 
and  contains  the  minimum  number  of  equations  to  compute 
a  given  set  of  outputs,  Yt.  As  a  result,  submodels  computed 
this  way  contain  only  a  small  subset  of  the  equations  of  the 
model  that  are  decoupled  from  the  rest  of  the  system  model 
equations.  In  general,  any  set  of  variables  in  the  system  can 
be  chosen  as  the  local  inputs  to  the  submodel,  Ui.  The  choice 
of  Ui  and  the  constraints  to  fulfill  depends  on  the  particular 
problem  to  be  solved  through  model  decomposition. 

The  first  model  decomposition  problem  computes  minimal 
submodels  from  the  nominal  system  model.  For  this  model 
decomposition  problem,  constraints  are  that  submodels  use 
the  global  model  inputs  and  some  measured  values  as  local 
inputs,  i.e.,  Ui  C  U  U  {Y  —  Yt).  An  algorithm  for  comput¬ 
ing  the  set  of  minimal  submodels  that  satisfies  these  prop¬ 
erties  is  given  in  (Daigle  et  ah,  2011),  which  is  based  on 
the  model  decomposition  algorithms  presented  in  (Pulido  & 
Alonso-Gonzalez,  2004;  Bregon  et  ah,  2012). 

The  second  model  decomposition  problem  computes  sub¬ 
models  for  residual  generation  and  fault  isolation.  For  this 
model  decomposition  problem,  the  constraints  are  that  sub¬ 
models  are  constructed  by  merging  the  minimal  submodels, 
to  fulfill  global  diagnosability  conditions  as  in  (Bregon  et 
ah,  2011). 3  Once  the  globally  diagnosable  subsystems  have 
been  designed,  the  merged  submodels  are  used  for  distributed 
residual  generation  and  to  compute  event-based  local  diag¬ 
noses  for  fault  isolation.  These  design  and  diagnoser  com¬ 
putation  processes  are  detailed  in  (Bregon  et  ah,  2011).  Be¬ 
cause  the  subsystems  are  designed  to  be  globally  diagnosable, 
the  resulting  local  diagnoses  are  independent,  and  can  pro¬ 
vide  globally  correct  diagnosis  results  without  a  centralized 
coordinator. 

The  third  model  decomposition  problem  computes  submod¬ 
els  for  distributed  fault  identification.  For  each  consistent 
fault  hypothesis  /,  the  joint  state -parameter  estimators  are 
computed  from  the  minimal  submodels  of  the  faulty  system 
model  M(f)  with  the  constraints  that  Ui  C  ULi(Y— Yj)  with 

-A  faulty  submodel  for  a  fault  /  is  denoted  as 

3  In  this  work,  a  subsystem  is  globally  diagnosable  if  all  faults  in  the  subsys¬ 
tem  are  distinguishable  from  every  other  fault  in  the  system  using  only  local 
measurements. 
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Distributed  Diagnosis  Distributed  Prognosis 


Figure  1.  An  instantiation  of  the  integrated  diagnosis  and  prognosis  architecture. 


Yt  as  a  singleton.  It  will  be  shown  later  that  the  fault  identi¬ 
fication  module  is  the  central  part  of  our  diagnosis-prognosis 
integration  approach  and  provides  the  joint  state-parameter 
estimations  for  the  prediction  module. 

Finally,  for  distributed  prediction,  the  model  decomposition 
problem  starts  off  from  the  faulty  system  model,  and,  as  de¬ 
tailed  in  (Daigle,  Bregon,  &  Roychoudhury,  2012),  it  fulfills 
the  following  constraints:  (i)  the  submodels  use  f/;  C  Up, 
where  Up  C  X  U  U  (here.  Up  is  a  set  of  variables  whose 
future  values  can  be  predicted  a  priori,  which  depends  on  the 
hypothesized  faults);  and  (ii)  each  computed  submodel  has 
at  least  one  c  £  C eo l$  belonging  to  C),  and  over  all  sub¬ 
models,  all  constraints  in  Ceols  are  covered.  This  ensures 
that  TpoLf  may  be  computed  for  the  system  from  the  local 
constraints. 

3.2.  Distributed  Architecture 

Figure  1  illustrates  an  example  architecture  for  our  distributed 
diagnosis  and  prognosis  scheme.  At  each  discrete  time  step, 
k,  the  system  takes  as  input  both  u/,  and  y/,:  and  splits  them 
into  local  inputs  uj.  and  local  outputs  y\.  for  the  local  diag¬ 
noses.  Within  each  AT;  local  diagnoser,  nominal  tracking  is 
performed,  computing  estimates  of  nominal  measurements, 
ylk.  The  fault  detector  compares  the  estimated  measurements 
against  the  observed  measurements,  to  determine  statistically 
significant  deviations  for  the  residual,  rk  =  yk  —  y* .  Qual¬ 
itative  values  of  the  deviations  in  the  residuals  are  used  by 
the  event-based  diagnoser  to  isolate  faults.  The  set  of  iso¬ 
lated  fault  candidates  Flk  together  with  the  estimated  nom- 

inal  states,  xj^,  parameters,  6k,  and  the  measurements,  yk , 
are  used  as  input  for  the  corresponding  identification  mod¬ 
ule.  Identification  is  performed  for  each  hypothesized  fault 
in  a  distributed  way,  e.g.,  for  the  isolated  faults  /i  and  f-2  in 
Figure  1,  we  run  an  instantiation  of  the  identification  sub¬ 
model  for  each  one  the  faults,  i.e.,  Xii(fi)  and  Xii(f-2)- 
Fault  identification  uses  the  minimal  submodels  from  the 
faulty  system  model,  and  computes  local  state-parameter  es¬ 


timates  p(xj  k,  Olf,k\ylo-k)-  These  local  estimates  are  then 
used  as  input  to  the  prediction  submodels.  In  some  cases, 
the  local  estimates  have  to  be  split  or  merged  with  other 
estimates  according  to  the  prediction  submodels.  For  ex¬ 
ample,  in  the  figure,  estimates  from  AIi(/i),  are  used  by 
both  local  prediction  submodels  Ad^/i)  and  Xi$(fi),  and 
those  submodels  may  also  need  estimates  not  included  within 
submodel  Adi(/i).  These  estimates  are  typically  obtained 
from  the  local  diagnosers  or  other  fault  identification  blocks.4 
Distributed  prediction  modules  compute,  for  each  hypothe¬ 
sized  fault,  local  EOL/RUL  predictions,  p^EOL^  kp  |y(j.;;.;, ) 
and  p{RUL)kp  |yj.fcp),  at  given  prediction  time  kp  based 
on  the  local  EOL  constraints.  Finally,  local  predictions 
are  combined  into  global  predictions  p(EOLjkp  |yo:fcj= )  and 
p(RULf  kp\y0.kp)  for  each  hypothesized  fault.  The  next 
section  describes  the  details  of  the  different  modules  of  the 
distributed  integrated  diagnosis  and  prognosis  architecture. 

4.  Diagnosis  and  Prognosis  Approach 

Figure  1  shows  the  basic  modules  of  our  distributed  integrated 
approach.  In  this  section  we  give  details  on  how  each  mod¬ 
ule  is  implemented,  and  establish  the  integration  between  the 
diagnosis  and  prognosis  tasks. 

4.1.  Distributed  Diagnosis 

For  distributed  diagnosis,  each  local  diagnoser  first  takes  a 
subset  of  the  local  inputs  uj,  and  local  outputs  ylk,  to  compute 
an  estimate  of  its  output  measurements  yk.  Tracking  is  per¬ 
formed  in  discrete  time  using  a  robust  filtering  scheme,  e.g., 
the  extended  or  unscented  Kalman  filter  (Julier  &  Uhlmann, 
2004),  which  provides  accurate  tracking  in  the  presence  of 
sensor  noise,  process  noise,  and  discretization  error. 


4Since  prediction  submodels  are  constructed  by  using  any  variable  which 
value  can  be  hypothesized  as  input,  in  some  cases,  prediction  submodels 
cannot  always  be  formed  by  just  merging  the  minimal  estimation  submod¬ 
els.  To  indicate  this,  we  named  our  prediction  submodels  differently  from 
the  estimation  submodels,  e.g.,  .M4 (/i )  and  .Ms (/i )  instead  of  A4 1  (/i )• 
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For  fault  detection,  a  statistical  test  is  used  to  look  for  signif¬ 
icant  deviations  in  the  residual  signal  rk,  which  is  computed 
as  the  difference  between  yj,  and  the  y‘k.  In  our  approach,  we 
use  a  Z-test  as  described  in  (Biswas  et  al.,  2003). 

Fault  isolation  is  performed  using  local  event-based  di¬ 
agnoses,  constructed  as  detailed  in  the  previous  sec¬ 
tion  (Bregon  et  ah,  2011;  Daigle  et  ah,  2009).  Fault  isolation 
is  triggered  when  a  fault  is  detected,  and  it  works  as  follows. 
Initially,  all  event-based  local  diagnoses  start  in  their  initial 
state,  and  the  set  of  faulty  candidates  is  empty.  Local  residual 
deviations  cause  the  local  diagnoses  to  move  from  one  state 
to  another.  These  residual  deviations  are  abstracted  to  a  tu¬ 
ple  of  qualitative  symbols  (<71,02)  for  each  residual  signal, 
where  a-\  represents  magnitude  changes  and  <72  represents 
slope  changes.  A  +  (resp.  — )  value  indicates  a  change  above 
(resp.  below)  normal  for  a  measurement  residual  or  a  positive 
(resp.  negative)  residual  slope.  A  0  implies  no  change  in  the 
measurement  value  or  a  flat  residual  slope.  The  symbols  are 
generated  using  a  sliding  window  technique  as  described  in 
detail  in  (Biswas  et  ah,  2003).  If  there  is  a  match  between 
an  event  from  the  current  state  and  a  tuple  of  qualitative  sym¬ 
bols  generated  by  any  residual,  the  local  diagnoser  moves  to 
the  next  state  and  remains  active.  If  not,  the  local  diagnoser 
blocks.  This  process  continues  until  a  local  diagnoser  reaches 
an  accepting  state,  which  corresponds  to  a  unique  isolation 
result. 

In  our  distributed  diagnosis  approach,  identification  submod¬ 
els,  A fj(/),  are  obtained,  as  explained  in  the  previous  sec¬ 
tion,  as  minimal  submodels  from  the  faulty  system  model.  A 
local  state -parameter  estimator  is  constructed  for  each  iden¬ 
tification  submodel  AL;(/),  and  produces  a  local  estimate 
p(x(-  k,  9 lf  k|yg.fc)  by  using  an  appropriate  algorithm.  In  this 
paper,  we  use  an  unscented  Kalman  filter  (UKF)  (Julier  & 
Uhlmann,  2004)  with  a  variance  control  algorithm  (Daigle, 
Saha,  &  Goebel,  2012). 

4.2.  Distributed  Prognosis 

The  local  state-parameter  estimates  for  each  local  distributed 
prediction  module  are  constructed  from  the  local  estimates  of 
the  distributed  fault  identification  submodels.  Each  predic¬ 
tion  submodel  is  made  up  of  a  set  of  states  and  parame¬ 
ters  ©j,  and  constructs  a  local  distribution p(x^  k,  9^  fc|yg.fc), 
by  assuming  that  the  local  state-parameter  estimates  are  suf¬ 
ficiently  represented  by  a  mean  /i'  and  covariance  £\  For 
each  prediction  submodel  Alj(/),  we  combine  the  estimates 
of  the  local  identification  submodels  that  estimate  states  and 
parameters  in  X,  U  0 ,  into  //'  and  S  '  .  If  two  submodels  esti¬ 
mate  the  same  state  variable  or  parameter,  then  many  different 
techniques  can  be  applied  depending  on  the  desired  perfor¬ 
mance  of  the  prediction  submodels,  e.g.,  taking  the  estimate 
with  the  smallest  variance,  or  taking  an  average. 

Several  approaches  can  be  used  to  perform  prediction  for  each 


Algorithm  1  EOL  Prediction 

Inputs:  {(x£\*g> 

Outputs:  {EOLk:^ , 

for  j  —  1  to  A  do 

k  <—  kp 

i(j)  iU) 

*-k  *-kp 

a‘U)  .  aiU) 

°k  bkp 

while  T'EOL(xlk), 

=  0do 

Predict  uj. 

~p(°i+ 1 

\oiU)) 

Xfc+1  ~  P(Xfc+ 1 

k  <—  k  +  1 

i(j)  i(j) 

k  -H+l 

piU) 
Wfc+ 1 

end  while 

EOLiu)  <-  k 

K  p 

end  for 

prediction  submodel.  In  this  work,  given  the  mean  and  co- 
variance  information,  we  represent  the  distribution  with  a  set 
of  sigma  points  derived  using  the  unscented  transform.  Then, 
each  sigma  point  is  simulated  forward  to  EOL,  and  we  re¬ 
cover  the  statistics  of  the  EOL  distribution  given  by  the  sigma 
points  (Daigle  &  Goebel,  2010). 

Algorithm  1  (Daigle,  Saha,  &  Goebel,  2012),  shows  the 
pseudocode  for  the  prediction  procedure.  The  algorithm 
is  executed  for  each  submodel  i,  deriving  local  EOL  pre¬ 
dictions  using  its  local  threshold  function  based  on  the  lo¬ 
cal  EOL  constraints.  For  a  given  submodel,  each  sample 
j  is  propagated  forward  until  T’EOL  (x 'j3k  ,  9 'jJk  )  evaluates 
to  1.  The  algorithm  hypothesizes  future  inputs  uk.  Then, 
the  global  EOL/RUL  is  determined  by  the  minimum  of  the 
local  EOL/RUL  distributions  for  each  prediction  submodel, 
i.e., p(EOLif  kp \yi0.kp )  andp(RULykp\ylkp).  Tocompute 
this,  we  sample  from  each  local  EOL  distribution  and  take 
the  minimum  of  the  local  samples.  This  is  repeated  many 
times  and  the  statistics  of  the  global  EOL  distribution  are 
computed  (Daigle,  Bregon,  &  Roychoudhury,  2012). 

5.  Case  Study 

In  this  section,  we  apply  our  distributed  diagnosis  and  prog¬ 
nosis  approach  to  a  four-wheeled  rover  testbed  developed  at 
NASA  Ames  Research  Center.  We  develop  a  model  of  the 
rover,  and  demonstrate  the  approach  using  simulated  scenar¬ 
ios. 

5.1.  Nominal  System  Modeling 

The  rover  model  was  originally  presented  in  (Balaban  et  al., 
2011).  In  this  section  we  summarize  the  main  features  and 
include  some  extensions  to  the  model. 

The  rover  consists  of  a  symmetric  rigid  frame  with  four 
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independently-driven  wheels.  The  wheel  speeds  are  governed 
by 

d>FL  =  -J- —  (rmFI  —  TfFL  -  rglFL  +TgrFL )  ,  (Cl) 

•Jfl 

d>FR  =  — -  —  TfFR  ~  TglFR  ~  TgrFR)  ,  fe) 

Jfr 

d 'BL  =  -7 {TmBL  ~  TfBL  ~  TglBL  +  TgrBL )  ,  (C3) 

■JBL 

^BR  =  {jmBR  —  TfBR  ~  TglFR  ~  TgrBR.)  •  (C4) 

Jbr 

The  F,  B,  L,  and  R  subscripts  stand  for  front,  left,  back, 
and  right,  respectively.  Here,  for  wheel  w,  ./,,  denotes  the 
wheel  inertia;  Tmw  is  the  motor  torque;  Tfw  =  Hfu,ojw  is 
the  wheel  friction  torque,  where  [i fw  is  a  friction  coefficient; 
Tgiw  =  Twiigi(vw  —  v)  is  the  torque  due  to  slippage,  where 
rw  is  the  wheel  radius,  figi  is  a  friction  coefficient,  vw  is  the 
translational  wheel  velocity,  and  v  is  the  translation  velocity 
of  the  rover  body;  and  Tgrw  =  rwngrwui  cos  7  is  the  torque 
due  to  the  rotational  movement  of  the  rover  body,  where  figrw 
is  a  friction  coefficient,  w  is  the  rotational  velocity  of  the  rover 
body,  and  7  =  arctan  l/b  with  l  being  the  rover’s  length  and 
b  being  its  width. 

The  translational  velocity  v  of  the  rover  is  described  by 

v  =  —  ( FglFL  +  FglFR  +  FglBL  +  FgiBR)  ,  (05) 

where  m  is  the  rover  mass,  and  for  wheel  w,  Fgiw  =  /igi  (vw  — 
v)  is  the  force  due  to  slippage.  The  rotational  velocity  w  is 
described  by 

UJ  =  j(d cos  jFgipR  +  d cos 'yFgisn  -  d cos ~iFglFL 

dcos'yFgiBp  dFgrFL  dFgrFR  dFgrBR 

-  dFgrBR).  (C6) 

Here,  J  is  the  rotational  inertia  of  the  rover  and  d  is  the  dis¬ 
tance  from  the  center  of  the  rover  to  each  wheel. 

The  wheels  are  driven  by  DC  motors  with  PI  control  that  sets 
the  voltages  V  applied  to  the  motors.  The  motor  currents  i 
are  governed  by 

IFL  =  -  ZFlRfL  -  kuU!FL),  (C7) 

ZFR  =  jr(VFR  -  1FrRFr  -  k^ulpn),  (C8) 

ZBL  =  rjr(VBL  ~  ZBlRbL  ~  kujlVBL),  (C9> 

i-BR  =  ~jr(VBR  -  IbrRbR  -  kujUJBR),  (cio) 

where  L  is  the  motor  inductance,  R  is  the  motor  resistance, 
and  ku  is  an  energy  transformation  term.  The  motor  torque 
is  Tmw  =  kTiw,  where  kT  is  an  energy  transformation  gain. 
The  voltages  applied  to  the  motors  are  determined  by  the  con¬ 
trollers,  where  for  wheel  w,Vw  =  P  *  (uw  —  ojw)  +  I  *  e,iw, 
where  P  is  a  proportional  gain,  uw  is  the  commanded  wheel 
speed,  I  is  an  integral  gain,  and  eiw  is  the  integral  error  term. 


The  integral  error  terms  are  governed  by 


eiFL  =  UFL  -  9JpB ,  (Cll) 

eiFR  =  uFR  -  u jfr,  (C12) 

eisz,  =  uBL  -  uBL,  (C13) 

eiBR  =  uBr  -  uiBR.  (C14) 


The  batteries,  which  are  connected  in  series,  are  described 
by  a  simple  electrical  circuit  equivalent  model  that  includes  a 
large  capacitance  Cb  in  parallel  with  a  resistance  Rp,  together 
in  series  with  another  resistance  Rs.5  The  battery  charge  vari¬ 


ables  qi  are  governed  by 

91  =  —Vi/Rpi  -  (iFL  +  iFR  +  iBR  +  zBl ),  (C15) 

92  =  —Vi/Rpi  —  (lFL  +  ZFR  +  ZBR  +  IBL),  (C16) 

93  =  -Vs/Rpz  -  ( ipL  +  iFR  +  iBR  +  ibl),  (C17) 

94  =  —Vi/Rpi  —  ( %FL  +  i-FR  +  ZBR  +  Zbl)-  (ci8) 


The  available  sensors  measure  the  voltages  of  the  batteries. 


Vi*  =  qi/Cbi  -  Rsi  *  ( iFL  +  zFR  +  zbr  +  zbl),  (C19) 

V2*  =  92/C62  -  Rs 2  *  (ZFL  +  ZFR  +  ZBR  +  Zbl),  (C20) 

V3*  =  93/C63  -  Rs 3  *  (zFL  +  zFR  +  iBR  +  zBl),  (C21) 

V4*  =  qi/Cbi  —  Rsi  *  (zfl  +  zfr  +  zBr  +  zbl),  (C22) 

the  motor  currents, 

Zfl  =  zfl,  (C23) 

z*fr  =  zfr,  (C24) 

zbl  =  zbl,  (C25) 

Z*BR  =  ZBR,  (C26) 

and  the  wheel  speeds, 

U>PL  =UFL,  (C27) 

U*PR  =  ^FR,  (C28) 

uBL=u}BL,  (C29) 

Z°*BR  =  UBR-  (C30) 

Here,  the  *  superscript  indicates  a  measured  value. 


5.2.  Faulty  System  Modeling 

In  this  work,  we  consider  different  faults  in  the  motors  and  the 
batteries.  First,  we  consider  friction-based  damage  progres¬ 
sion  in  the  motors,  resulting  in  an  increase  in  motor  friction 
over  time.  For  wheel  w,  the  fault  progression  function  is  de¬ 
fined  as: 

(zfFL  =  ufFL  fZfFL  UPL’  (c3l) 

hfFR  =  FfpR  RfPR  UpR,  (C32) 

hfBL  =  VfBL  tZfBL  UBL,  (c33) 

t ZfBR  =  VfBR  RfBR  UBR,  (C34) 


5  We  use  a  simple  model  here  only  for  demonstration  purposes.  More  detailed 
battery  models  for  prognosis  can  be  found  in  the  literature,  e.g.,  (Saha  & 
Goebel,  2009). 
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Submodel 

Xi 

<0/ 

M7i 

Yi 

<Ji 

Mi 

91 

Cbi 

lFL ,  lFR,  %BL,%BR 

Cl5,Cl9,C23,C24,C25,C26 

M2 

92 

cb2 

%FL ,  lFR,  lBL,lBR 

v2* 

Cl6,C20,C23,C24,C25,C26 

Ms 

93 

Cb3 

I'FL ,  lFR,  %BL,%BR 

V3* 

Cl7,C21,C23,C24,C25,C26 

Ma 

94 

Cbi 

%FL ,  lFR,  iBL,1’BR 

VC 

Cl8,C22,C23,C24,C25,C26 

Ms 

IFL,  e-iFL 

Rfl 

UfL,U1*fl 

i*FL 

C7,Cn,C23,C27 

Ms 

IFR,  e-iFR 

Rfr 

ufr,  ui’fr 

i*FR 

C8,Cl2,C24,C28 

M-7 

iBL,e.iBL 

Rbl 

UbL,U1*bl 

i*BL 

C9,Cl3,C25,C29 

Ms 

IBR ,  e-iBR 

Rbr 

UBR,  UIrr 

Cl0,Cl4,C26,C30 

Mg 

UJFL,V,  Ul,  [IfFL 

VfFL 

'I’FL,UfR,UbL,UbR 

U*FL 

Cl  ,C31  ,C5  ,C6,C23  ,C28,C29  ,C30 

M 10 

U)FR,  V,  Ul,  fifFR 

RfFR 

\fr, ufl , uBl , uBr 

Ufr 

C2,C32,C5,C6,C24,C27,C29,C30 

Atn 

U>BL,V,  Ul,  fifBL 

VfBL 

lBL,UpL,UFR ,  UrR 

u*bl 

C3,C33,C5,C6,C25,C27,C28,C30 

At  12 

U>BR,  V,  Ul,  fifBR 

VfBR 

^BR,Ufl,UFR,UBL 

Ubr 

C4,C34,C5,C6,C26,C27,C28,C29 

Table  1 .  Set  of  minimal  submodels  for  the  rover  testbed  computed  from  the  nominal  system  model. 

Submodel 

Aj 

Ui 

Yi  Ui 

Mr,, 9 

c VFL,V,UJ, 

[IfFL ,  IFL,  eiFL 

VfFL,RFL  UfL,U>fr,U1bl,Ubr 

u’fl  ,  iFL  Us  U  U9 

Ate, 10 

uiFr,v,uj, 

flfFR,  IFR,  e.iFR 

VfFR,RFR  UfR,U*Fl,U*bl,U*bR 

Ufr,  ^fr  Cs  U  C10 

At7.i1 

ujbl,v,uj, 

[AfBL ,  iBL,  e,BL 

RfBL,  Rbl  ubl,ui*fl,uj*fr,lo*Br 

U*BL,i*BL  C7  U  Cll 

At8,12 

lobr,  v,  Ul, 

flfBR,  IBR,  e.iBR 

RfBR,  Rbr  ubr,Ufl,Ufr,Ubl 

Ubr,Ybr  Us  U  C12 

At  1,2, 3, 4 

91,92,93,94 

Cbl,Cb2,Cb3,Cb4  IflYfrYblUbR 

Vi*,  V2*,  V3*,  VC  C1UC2UC3U  U4 

Table  2.  Residual  generation  and  fault  isolation  submodels. 


where  /ifw  is  the  fault  parameter,  and  Ufw  is  the  fault  pro¬ 
gression  parameter. 

We  also  consider  abrupt  resistance  increases  in  the  motors, 
represented  as  an  abrupt  change  in  parameter  Rw  for  wheel 
to,  with  A Rw  as  the  fault  progression  parameter. 

For  the  batteries,  we  consider  abrupt  capacitance  decreases, 
represented  as  an  abrupt  change  in  parameter  6'/,,  for  capacity 
i.  A Cbi  is  the  fault  progression  parameter. 

We  are  interested  in  predicting  when  any  of  the  rover  batteries 
are  at  their  charge  threshold,  beyond  which  the  batteries  will 
be  damaged.  These  faults  can  cause  the  charge  thresholds  to 
be  reached  earlier  since  they  will  affect  current  draw.  The 
constraints  are  given  as 

91  >  9“ ,  P-35 ) 

92  >  q~ ,  (C36) 

93  >  q~  ,  (C37) 

94  >  q~  ,  (C38) 

where  the  charge  threshold  is  given  by  q~  =  2  x  104  C.  The 
rover  cannot  be  operated  when  any  of  the  constraints  C35-C38 
are  violated. 

5.3.  Results 

To  demonstrate  the  validity  of  the  approach,  we  describe  two 
different  faulty  scenarios  of  the  rover.  In  the  first,  friction 
damage  is  progressing  on  one  motor,  and  in  the  second,  a 
capacitance  decrease  occurs  in  one  battery.  In  all  cases,  the 
rover  travels  between  various  waypoints,  moving  at  an  aver¬ 
age  speed  of  0.5  m/s.  Table  1  shows  the  minimal  submodels 
for  the  rover  derived  by  using  measured  values  as  local  in¬ 
puts.  Table  2  shows  the  submodels  for  residual  generation 


Submodel 

Xi 

AT 

Ui 

Y 

At  17(0,1 ) 

9i,  Cbi 

“aot 

I’FL,  ^FR,  %BL,iBR 

0 

Cl5,Cl9,C35 

Mis{Ub2) 

92,  Cb  2 

AC„2 

iFL,lFR,  %BL,iBR 

0 

Cl6,C20,C36 

At  19  (Os) 

93,  Cb 3 

AC63 

iFL,lFR,  %BL,iBR 

0 

Cl7,C21,C37 

At  20(0,4) 

94,  Cm 

ACm 

lFL,iFR,  %BL,iBR 

0 

Cl8,C22,C38 

Table  4.  Prediction  submodels  for  capacitance  faults. 


and  fault  isolation.  These  submodels  have  been  designed  to 
obtain  globally  diagnosable  subsystems  by  using  the  design 
algorithm  in  (Bregon  et  al.,  201 1).  In  this  work,  we  have  con¬ 
sidered  five  subsystems,  one  for  each  wheel  components  and 
another  one  for  the  batteries.  For  example,  the  subsystem  for 
the  front  left  wheel  components  is  not  globally  diagnosable  if 
we  only  consider  submodel  Atg  (which  includes  the  front  left 
wheel  friction  wear  parameter,  v/fl)-  The  design  algorithm 
in  (Bregon  et  al.,  2011)  determines  that  we  need  to  merge 
submodels  Ad  5  and  Ad  9  to  make  the  front  left  wheel  subsys¬ 
tem  globally  diagnosable.  The  process  is  similar  for  the  rest 
of  the  subsystems. 

Minimal  submodels  for  identification  can  be  computed  from 
the  minimal  submodels  in  Table  1  by  defining  the  fault  pro¬ 
gression  function  (if  necessary),  and  by  making  the  fault  pa¬ 
rameter  to  become  a  state  and  the  corresponding  fault  pro¬ 
gression  parameter  to  become  the  parameter.  Regarding  pre¬ 
diction,  the  correct  prediction  submodels  to  use  depend  on  the 
scenario,  as  will  be  shown  later. 

5.3.1.  Friction  Damage  Progression 

We  first  consider  a  scenario  in  which,  for  the  front-left  mo¬ 
tor,  the  friction  begins  to  increase.  The  friction  damage  pro¬ 
gression  begins  at  t  =  50  s  with  friction  wear  parameter 
v/FL  =  1  x  10-3  s.  A  fault  is  detected  by  the  local  di- 
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Submodel  Xj  0)  Ul  Y~Oi 

■Ml3([AfFL)  qi,iFL,  dFL ,  Wfi ,  HfFL  Cfcl,  vfFL,  RfL  UFL ,  V ,  U} ,  iFR,  iBL ,  iBR  0  Cl  ,C31  ,C7,Cll  ,Cl5  ,C35 

■Mlii/AfFL)  q2,  iFLi  ZiFL,  UfL,  VfFL  Cb2,  VfFL ,  RfL  UfL,V  ,  (J  ,iFR,iBL,iBR  0  Cl,C3i,C7,Cn,Cl6,C36 

MlbiUfFL.)  <?3,  IFL,  eiFL,  UFL,  HfFL  Cb3 ,  V  f  F  L ,  RfL  UfL  ,  V ,  LJ ,  iFR,  iBL ,  iBR  0  Cl,C31,C7,Cll,Cl7,C37 

M 16  (/l  fFL)  <?4,  IFL,  eiFL,  UIFL,  ^fFL  CbH,  Vj  F  L ,  RfL  UFL  ,  V ,  UJ ,  tFfl,  iBL ,  is  R  0  Cl,C3i,C7,Cll,Cl8,C38 


Table  3.  Prediction  submodels  using  commanded  wheel  speeds  and  rover  velocities  as  local  inputs. 


Figure  2.  Estimated  VfFL  values. 
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Figure  3.  Current  i*FL  increase  through  time. 


agnoser  computed  from  submodel  A4r,s>  at  119.25  s,  via  an 
increase  in  the  motor  current  ipL-  The  initial  candidate  list  is 
immediately  reduced  to  one  candidate,  {zz/fl},  based  on  the 
signatures  and  orderings  (other  faults  in  the  front  left  wheel, 
like  Rfl,  produce  different  fault  signatures).  Thus  the  true 
fault  is  isolated. 


incorrect  to  use  it  as  a  local  input  for  prediction  and  to  de¬ 
compose  the  prediction  problem  into  independent  local  pre¬ 
diction  problems  for  the  batteries  and  motors,  i.e.,  it  is  not 
known  a  priori.  Therefore,  we  compute  submodels  using  as 
local  inputs  average  values  for  the  remaining  motor  currents, 
average  commanded  wheel  speeds,  and  average  rover  trans¬ 
lational  velocity  v  and  rotational  velocity  u>.  The  prediction 
submodels  for  this  case  are  shown  in  Table  3.  EOL  for  this 
fault  is  computed  by  merging  the  local  EOL  from  those  sub¬ 
models  in  the  table.  Note  that  the  prediction  submodels  used 
in  this  case  do  not  correspond  directly  to  those  used  for  es¬ 
timation.  So,  when  constructing  the  estimate  for  Ad  13,  for 
example,  it  takes  the  estimates  from  Adi  and  Ad 9. 

The  prediction  results  are  shown  in  Figure  4.  The  increased 
friction  causes  the  batteries  to  discharge  faster,  and  EOL  oc¬ 
curs  around  1650  s.  Here,  we  used  the  relative  accuracy  (RA) 
as  a  measure  of  prediction  accuracy,  and  the  relative  standard 
deviation  (RSD)  as  a  measure  of  spread.  Each  prediction  met¬ 
ric  is  averaged  over  multiple  prediction  points  (one  every  100 
s  of  usage)  (see  (Saxena  et  al.,  2010;  Daigle,  Saha,  &  Goebel, 
2012)  for  the  mathematical  definitions  of  these  metrics).  For 
this  experiment,  RA  averages  to  91.63%  and  RSD  averages 
to  16.26%. 

For  the  sake  of  comparison,  we  also  ran  this  experiment  us¬ 
ing  the  centralized  approach.  Figure  5  shows  the  prediction 
results  obtained.  Looking  at  the  prediction  metrics,  we  see 
that  the  centralized  approach  behaved  very  similar  to  the  dis¬ 
tributed  approach  but  a  little  bit  worse,  with  RA  averaging 
90.90%  and  RSD  averaging  17.72%.  However,  this  is  just 
a  particular  example,  but,  in  general,  both  approaches  obtain 
equivalent  results. 

5.3.2.  Capacitance  Decrease 


Fault  identification  is  initiated  once  the  candidate  is  isolated. 
For  the  friction  damage  progression  fault,  the  wear  rate  zz/fl 
estimate  averages  to  v/fl  =  1  x  10~3  s  with  very  small 
output  error.  Figure  2  shows  the  wear  parameter  estimate  for 
friction  damage. 

As  a  result  of  the  continuously  increasing  friction,  the  current 
drawn  by  the  motor  increases  as  well  in  order  for  the  motor 
controller  to  maintain  the  same  desired  wheel  speed  (Figure 
3  shows  this  increase  in  the  current  through  time).  Hence,  the 
total  current  drawn  from  the  batteries  is  increased,  and  EOL 
occurs  around  half  an  hour.  Because  i fl  is  constantly  chang¬ 
ing,  and  in  a  way  that  is  dependent  on  the  motor  state,  it  is 


As  a  second  scenario,  we  consider  a  capacitance  decrease 
fault  in  battery  3  of  the  rover,  Cb3 .  The  fault  begins  at  t  =  50 
s  with  an  abrupt  decrease  from  2000  to  1800  in  the  capacity 
of  the  battery.  The  fault  is  detected  immediately  by  the  local 
diagnoser  computed  from  submodel  Ad  1,2,3, 4  at  50.0  s,  via  an 
increase  in  the  voltage  V3.  The  fault  candidate  is  immediately 
isolated,  { Ct,:', } ,  based  on  the  signatures  and  orderings,  thus 
starting  the  fault  identification.  For  the  capacitance  fault,  the 
estimated  value  of  the  capacitance  averaged  Cb3  =  1798.6 
C  with  very  small  output  error.  As  a  result  of  the  decrease 
in  capacitance,  the  battery  discharges  at  a  faster  rate,  and  so 
reaches  end  of  discharge  more  quickly.  The  prediction  sub- 
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Figure  4.  Predicted  RUL  of  the  rover  for  the  distributed  ap¬ 
proach.  The  mean  is  indicated  with  a  dot  and  confidence  in¬ 
tervals  for  5%  and  95%  by  lines.  The  gray  cone  depicts  an 
accuracy  requirement  of  15%. 
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Figure  5.  Predicted  RUL  of  the  rover  for  the  centralized  ap¬ 
proach.  The  mean  is  indicated  with  a  dot  and  confidence  in¬ 
tervals  for  5%  and  95%  by  lines.  The  gray  cone  depicts  an 
accuracy  requirement  of  15%. 


models  for  faults  in  the  capacity  of  the  batteries  are  shown 
in  Table  4.  For  this  scenario,  with  a  fault  in  Cb3 ,  we  used 
submodel  yVLig,  obtaining  RA  average  to  98.25%  and  RSD 
average  to  10.12%. 

6.  Conclusions 

This  paper  presented  a  distributed  integrated  model-based 
diagnosis  and  prognosis  framework.  Our  approach  starts 
off  with  a  common  modeling  paradigm  to  model  both  the 
nominal  behavior  and  fault  progression,  and  then  proposes 
a  framework  where  the  global  system  model  is  decomposed 
into  smaller  independent  submodels.  These  submodels  are 
then  used  to  distribute  the  different  diagnosis  and  prognosis 
tasks.  Model  decompostion  is  carried  out  based  on  the  re¬ 
quirements  and  constraints  of  each  task.  We  demonstrated 
our  approach  on  a  four-wheeled  rover  testbed,  where  we  di¬ 
agnosed  faults  and  prognosed  the  EOL/RUL  accurately.  We 
compared  results  obtained  by  using  our  distributed  approach 
against  those  obtained  using  a  centralized  approach,  showing 
that  both  approaches  obtain  the  same  results. 

Most  approaches  in  the  literature  focus  in  either  the  diagnosis 
or  the  prognosis  task.  Some  works  have  proposed  the  inte¬ 
gration  of  both  tasks  within  a  common  framework  (Patrick 
et  al.,  2007;  Orchard  &  Vachtsevanos,  2009;  Roychoudhury 
&  Daigle,  2011),  however,  unlike  our  approach,  these  ap¬ 
proaches  perform  the  diagnosis  and  prognosis  tasks  in  a  cen¬ 
tralized  way,  thus  suffering  from  scalability  issues  due  to  the 
large  number  of  states  and  parameters  in  real-world  systems. 
To  the  best  of  our  knowledge,  there  is  no  approach  in  the  lit¬ 
erature  which  combines,  in  a  distributed  way,  the  diagnosis 
and  prognosis  tasks.  Our  approach  is  limited  by  the  number 
and  location  of  the  sensors  in  the  system.  Since  our  decom¬ 
position  algorithm  is  guided  by  the  set  of  available  sensors, 
the  distribution  capabilities  of  the  approach  is  determined  by 
them. 

In  future,  we  will  apply  this  approach  to  larger  systems,  to 
study  the  scalability  of  our  diagnosis  and  prognosis  scheme, 
and  will  perform  a  more  detailed  comparison  against  the  re¬ 
sults  obtained  by  using  a  centralized  approach.  We  will  also 
extend  the  capability  of  this  approach  to  hybrid  systems,  as 
well  as  diagnosis  and  prognosis  of  multiple  faults. 

Acknowledgments 

A.  Bregon’s  funding  for  this  work  was  provided  by  the  Span¬ 
ish  MCI  TIN2009- 11326  grant.  M.  Daigle  and  I.  Roychoud- 
hury’s  funding  for  this  work  was  provided  by  the  NASA 
System-wide  Safety  and  Assurance  Technologies  (SSAT) 
Project. 


9 


Annual  Conference  of  the  Prognostics  and  Health  Management  Society  2012 


References 

Balaban,  E.,  Narasimhan,  S.,  Daigle,  M.,  Celaya,  J.,  Roy- 
choudhury,  I.,  Saha,  B.,  et  al.  (2011,  September). 
A  Mobile  Robot  Testbed  for  Prognostics-Enabled  Au¬ 
tonomous  Decision  Making.  In  Annual  Conference 
of  the  Prognostics  and  Health  Management  Society 
(p.  15-30).  Montreal,  Canada. 

Biswas,  G.,  Simon,  G.,  Mahadevan,  N.,  Narasimhan,  S., 
Ramirez,  J.,  &  Karsai,  G.  (2003,  June).  A  robust 
method  for  hybrid  diagnosis  of  complex  systems.  In 
Proceedings  of  the  5th  Symposium  on  Fault  Detection, 
Supervision  and  Safety  for  Technical  Processes  (pp. 
1125-1131). 

Blanke,  M.,  Kinnaert,  M.,  Lunze,  J.,  &  Staroswiecki, 
M.  (2006).  Diagnosis  and  Fault-Tolerant  Control. 
Springer. 

Bregon,  A.,  Biswas,  G.,  &  Pulido,  B.  (2012).  A  Decompo¬ 
sition  Method  for  Nonlinear  Parameter  Estimation  in 
TRANSCEND.  IEEE  Trans.  Syst.  Man.  Cy.  Part  A, 
42(3),  751-763. 

Bregon,  A.,  Daigle,  M.,  Roychoudhury,  I„  Biswas,  G.,  Kout- 
soukos,  X.,  &  Pulido,  B.  (2011,  Oct).  Improving 
Distributed  Diagnosis  Through  Structural  Model  De¬ 
composition.  In  Proceedings  of  the  22nd  International 
Workshop  on  Principles  of  Diagnosis  (p.  195-202). 
Murnau,  Germany. 

Daigle,  M.,  Bregon,  A.,  &  Roychoudhury,  I.  (2011,  Septem¬ 
ber).  Distributed  Damage  Estimation  for  Prognostics 
Based  on  Structural  Model  Decomposition.  In  Pro¬ 
ceedings  of  the  Annual  Conference  of  the  Prognostics 
and  Health  Management  Society  201 1  (p.  198-208). 

Daigle,  M.,  Bregon,  A.,  &  Roychoudhury,  I.  (2012).  Dis¬ 
tributed  prognostics  based  on  structural  model  decom¬ 
position.  (Manuscript  submitted  for  publication) 

Daigle,  M.,  &  Goebel,  K.  (2010,  October).  Improv¬ 
ing  Computational  Efficiency  of  Prediction  in  Model- 
based  Prognostics  Using  the  Unscented  Transform.  In 
Annual  Conf.  of  the  Prognostics  and  Health  Manage¬ 
ment  Society  2010. 

Daigle,  M„  Koutsoukos,  X.,  &  Biswas,  G.  (2009,  July).  A 
Qualitative  Event-based  Approach  to  Continuous  Sys¬ 
tems  Diagnosis.  IEEE  Trans,  on  Control  Systems  Tech¬ 
nology,  77(4),  780-793. 

Daigle,  M.,  Saha,  B„  &  Goebel,  K.  (2012,  March). 
A  comparison  of  filter-based  approaches  for  model- 
based  prognostics.  In  Proceedings  of  the  2012  IEEE 
Aerospace  Conference. 

Gertler,  J.  J.  (1998).  Fault  Detection  and  Diagnosis  in  Engi¬ 
neering  Systems.  New  York,  NY:  Marcel  Dekker,  Inc. 

Julier,  S.  J.,  &  Uhlmann,  J.  K.  (2004,  March).  Unscented 
filtering  and  nonlinear  estimation.  Proc.  of  the  IEEE, 
92(3),  401^122. 

Luo,  J.,  Pattipati,  K.  R.,  Qiao,  L.,  &  Chigusa,  S.  (2008, 


September).  Model-based  prognostic  techniques  ap¬ 
plied  to  a  suspension  system.  IEEE  Trans,  on  Systems, 
Man  and  Cybernetics,  Part  A:  Systems  and  Humans, 
38(5),  1156-1168. 

Mosterman,  P.  J.,  &  Biswas,  G.  (1999,  November).  Diag¬ 
nosis  of  Continuous  Valued  Systems  in  Transient  Op¬ 
erating  Regions.  Systems,  Man  and  Cybernetics,  Part 
A:  Systems  and  Humans,  IEEE  Trans,  on,  29(6),  554  - 
565. 

Orchard,  M.  E.,  &  Vachtsevanos,  G.  (2009).  A  Particle- 
filtering  Approach  for  On-line  Fault  Diagnosis  and 
Failure  Prognosis.  Trans,  of  the  Institute  of  Measure¬ 
ment  and  Control,  57(3/4),  221-246. 

Patrick,  R.,  Orchard,  M.  E„  Zhang,  B.,  Koelemay,  M., 
Kacprzynski,  G.,  Ferri,  A.,  et  al.  (2007,  September). 
An  Integrated  Approach  to  Helicopter  Planetary  Gear 
Fault  Diagnosis  and  Failure  Prognosis.  In  Proc.  of  the 
42nd  Annual  Systems  Readiness  Technology  Conf.  Bal¬ 
timore,  MD,  USA. 

Pulido,  B.,  &  Alonso-Gonzalez,  C.  (2004).  Possible  Con¬ 
flicts:  a  compilation  technique  for  consistency-based 
diagnosis.  IEEE  Trans,  on  Systems,  Man,  and  Cyber¬ 
netics,  Part  B,  Special  Issue  on  Diagnosis  of  Complex 
Systems,  34(5),  2192-2206. 

Roychoudhury,  I.  (2009).  Distributed  Diagnosis  of  Continu¬ 
ous  Systems:  Global  Diagnosis  Through  Local  Analy¬ 
sis.  Unpublished  doctoral  dissertation,  Vanderbilt  Uni¬ 
versity. 

Roychoudhury,  I.,  &  Daigle,  M.  (2011,  October).  An  inte¬ 
grated  model-based  diagnostic  and  prognostic  frame¬ 
work.  In  Proceedings  of  the  22nd  International  Work¬ 
shop  on  Principles  of  Diagnosis  (p.  44-51). 

Saha,  B.,  &  Goebel,  K.  (2009,  September).  Modeling  Li-ion 
battery  capacity  depletion  in  a  particle  filtering  frame¬ 
work.  In  Proc.  of  the  Annual  Conf.  of  the  Prognostics 
and  Health  Management  Society  2009. 

Saxena,  A.,  Celaya,  J.,  Saha,  B.,  Saha,  S.,  &  Goebel,  K. 
(2010).  Metrics  for  offline  evaluation  of  prognostic 
performance,  hit.  Journal  of  Prognostics  and  Health 
Management. 

Staroswiecki,  M.,  &  Declerck,  P.  (1989,  July).  Analyti¬ 
cal  redundancy  in  nonlinear  interconnected  systems  by 
means  of  structural  analysis.  In  I  FA  C  Symp.  on  Ad¬ 
vanced  Information  Processing  in  Automatic  Control. 

Williams,  B.,  &  Millar,  B.  (1998).  Decompositional  Model- 
based  Learning  and  its  analogy  to  diagnosis.  In  Proc.  of 
the  Fifteenth  National  Conf.  on  Artificial  Intelligence, 
AAAF98  (p.  197-204). 

Biographies 

Anibal  Bregon  received  his  B.S.,  M.S.,  and  Ph.D.  degrees  in 

Computer  Science  from  the  University  of  Valladolid,  Spain, 

in  2005,  2007,  and  2010,  respectively.  From  September 

2005  to  June  2010,  he  was  Graduate  Research  Assistant 


10 


Annual  Conference  of  the  Prognostics  and  Health  Management  Society  2012 


with  the  Intelligent  Systems  Group  at  the  University  of  Val¬ 
ladolid,  Spain.  He  has  been  visiting  researcher  at  the  Insti¬ 
tute  for  Software  Integrated  Systems,  Vanderbilt  University, 
Nashville,  TN,  USA;  the  Dept,  of  Electrical  Engineering, 
Linkoping  University,  Linkoping,  Sweden;  and  the  Diagnos¬ 
tics  and  Prognostics  Group,  NASA  Ames  Research  Center, 
Mountain  View,  CA,  USA.  Since  September  2010,  he  has 
been  Assistant  Professor  and  Research  Scientist  at  the  De¬ 
partment  of  Computer  Science  from  the  University  of  Val¬ 
ladolid. 

Dr.  Bregon  is  a  member  of  the  Prognostics  and  Health  Man¬ 
agement  Society  and  the  IEEE.  His  current  research  inter¬ 
ests  include  model-based  reasoning  for  diagnosis,  prognos¬ 
tics,  health-management,  and  distributed  diagnosis  of  com¬ 
plex  physical  systems. 

Matthew  Daigle  received  the  B.S.  degree  in  Computer  Sci¬ 
ence  and  Computer  and  Systems  Engineering  from  Rensse¬ 
laer  Polytechnic  Institute,  Troy,  NY,  in  2004,  and  the  M.S. 
and  Ph.D.  degrees  in  Computer  Science  from  Vanderbilt  Uni¬ 
versity,  Nashville,  TN,  in  2006  and  2008,  respectively. 

From  September  2004  to  May  2008,  he  was  a  Graduate  Re¬ 
search  Assistant  with  the  Institute  for  Software  Integrated 
Systems  and  Department  of  Electrical  Engineering  and  Com¬ 
puter  Science,  Vanderbilt  University,  Nashville,  TN.  During 
the  summers  of  2006  and  2007,  he  was  an  intern  with  Mission 
Critical  Technologies,  Inc.,  at  NASA  Ames  Research  Center. 
From  June  2008  to  December  201 1,  he  was  an  Associate  Sci¬ 


entist  with  the  University  of  California,  Santa  Cruz,  at  NASA 
Ames  Research  Center.  Since  January  2012,  he  has  been  with 
NASA  Ames  Research  Center  as  a  Research  Computer  Sci¬ 
entist.  His  current  research  interests  include  physics-based 
modeling,  model-based  diagnosis  and  prognosis,  simulation, 
and  hybrid  systems. 

Dr.  Daigle  is  a  member  of  the  Prognostics  and  Health  Man¬ 
agement  Society  and  the  IEEE.  He  is  a  recipient  of  a  Univer¬ 
sity  Graduate  Fellowship  from  Vanderbilt  University,  a  best 
paper  award  in  the  Annual  Conference  of  the  Prognostics 
and  Health  Management  Society  2011,  a  NASA  Ames  Group 
Achievement  Award  in  2011,  and  an  Ames  Contractor  Coun¬ 
cil  Excellence  Award  in  201 1 .  He  has  published  over  40  peer- 
reviewed  papers  in  the  area  of  Systems  Health  Management. 

Indranil  Roychoudhury  received  the  B.E.  (Hons.)  degree 
in  Electrical  and  Electronics  Engineering  from  Birla  Institute 
of  Technology  and  Science,  Pilani,  Rajasthan,  India  in  2004, 
and  the  M.S.  and  Ph.D.  degrees  in  Computer  Science  from 
Vanderbilt  University,  Nashville,  Tennessee,  USA,  in  2006 
and  2009,  respectively.  Since  August  2009,  he  has  been  with 
SGT,  Inc.,  at  NASA  Ames  Research  Center  as  a  Computer 
Scientist. 

Dr.  Roychoudhury  is  a  member  of  the  Prognostics  and  Health 
Management  Society  and  the  IEEE.  His  research  interests 
include  hybrid  systems  modeling,  model-based  diagnostics 
and  prognostics,  distributed  diagnostics  and  prognostics,  and 
Bayesian  diagnostics  of  complex  physical  systems. 


11 


